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Onconase® (ONC) is an amphibian member of the pancreatic ribonuclease 
superfamily that is selectively toxic to tumor cells. It is a much less efficient 
enzyme than the archetypal ribonuclease A and, in an attempt to gain fur- 
ther insight, we report the first atomic resolution crystal structure of ONC, 
determined in complex with sulfate ions at 100 K. The electron density 
map is of a quality sufficient to reveal significant nonplanarity in several 
peptide bonds. The majority of active site residues are very well defined, 
with the exceptions being Lys31 from the catalytic triad and Lys33 from 
the Bi subsite, which are relatively mobile but rigidify upon nucleotide 
binding. Cryocooling causes a compaction of the unit cell and the protein 
contained within. This is principally the result of an inward movement of 
one of the lobes of the enzyme (lobe 2), which also narrows the active site 
cleft. Binding a nucleotide in place of sulfate is associated with an approxi- 
mately perpendicular movement of lobe 2 and has little further effect on 
the cleft width. Aspects of this deformation are present in the principal 
axes of anisotropy extracted from C™ atomic displacement parameters, indi- 
cating its intrinsic nature. The three lowest-frequency modes of ONC 
motion predicted by an anisotropic network model are compaction/expan- 
sion variations in which lobe 2 is the prime mover. Two of these have high 
similarity to the cryocooling response and imply that the essential 'breath- 
ing' motion of ribonuclease A is conserved in ONC. Instead, shifts in con- 
formational equilibria may contribute to the reduced ribonucleolytic 
activity of ONC. 



Database 

Structural data have been submitted to the Protein Data Bank under accession number 3SNF. 



Introduction 

Onconase® (ONC, Alfacell Corporation, Somerset, NJ, 
USA) is a ribonuclease with anti-tumor activity iso- 
lated from oocytes of the Northern leopard frog, 
Rana pipiens [1]. It is selectively toxic toward numer- 
ous cancer cell types [2] and represents a promising 
anti-cancer drug, having undergone a confirmatory 



Phase Illb clinical trial for the treatment of nonresec- 
table malignant mesothelioma [3]. 

ONC is a member of the pancreatic ribonuclease 
superfamily whose archetypal member, ribonuclease A 
(RNase A; EC 3.1.27.5 ), holds a central position in the 
study of protein chemistry, enzymology and molecular 



Abbreviations 

ADP, atomic displacement parameter; ANM, anisotropic network model; ONC, Onconase; PDB, Protein Data Bank; RNase A, ribouclease A. 
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evolution [4-6]. Although ONC has only 30% 
sequence identity with RNase A, contains 20 fewer res- 
idues (104 versus 124) and possesses an alternative 
disulfide bonding arrangement, it has a highly similar 
cx/p fold [7], The structure is bi-lobed, with the lobes 
founded on P-sheets designated Vi and V 2 , respectively 
[8,9] (Fig. 1). Canonically, the active site lies in the 
cleft between the two lobes [10]. 

The ribonucleolytic activity of ONC is essential to 
its cytotoxicity, and the target of this activity is intra- 
cellular RNA [11]. In common with RNase A, RNA 
cleavage is endonucleolytic, occuring specifically on the 
3'-phosphate side of pyrimidine residues [12,13]. How- 
ever, whereas RNase A is an efficient enzyme exhibit- 
ing k cat /K m values in excess of 10 7 m _1 -s _1 with 
optimized oligonucleotide substrates [14], ONC is 
much less so under equivalent conditions (by ~ 10 3 - 
fold) [13]. Research to date suggests that the physical 
basis for this derives not only from localized structural 
features, but also from differences in the collective 
dynamics of the two enzymes. The body of work on 
RNase A is extensive [6] and has shown that RNase A 
interacts with the central dinucleotide of the substrate 
via multiple interactions spread across a tightly-fitting 
binding site [15,16]. Cleavage occurs across the P-O 5 
bond by a two-step general acid-base mechanism 
involving a catalytic triad in which two histidine resi- 




Fig. 1. ONC topology. Schematic representation of ONC bound to 
a dinucleotide substrate analogue as derived from the crystal struc- 
ture of the ONC-d(AUGA) complex (PDB entry 2I5S ) [10]. Helices 
H1-H3 and p-strands B0-B7 are labeled, and are numbered for con- 
sistency with the structural elements of RNase A (a short loop is 
found in place of p-strand B2 of RNase A). The foundations of 
lobes 1 and 2 are p-sheets V-, (B3, B6 and B7 comprising residues 
55-58, 86-91 and 96-101, respectively) and V 2 (B0, B1, B4 and B5 
comprising residues 11-12, 33-38, 63-70 and 77-84, respectively) 
[8,9]. The active site lies in a cleft between the lobes, here occu- 
pied by a d(UpGp) moiety shown in stick representation. The four 
disulfide bonds are shown in ball-and-stick representation (the one 
connecting H3 and B6 is obscured). 



dues (His 12 and His 119) orchestrate the movement of 
protons [17], and a lysine residue (Lys41) stabilizes the 
transition state via donation of a hydrogen bond to a 
nonbridging phosphoryl oxygen [18]. In addition, 
RNase A is known to interconvert between open and 
closed forms [8,19-24]. The motion is subtle but 
appears to be coupled to turnover [25], with an open 
conformation associated with substrate encounter and 
a closed one with catalysis [24,26]. Data relating to 
ONC are less abundant, although the enzyme is known 
to differ from RNase A in several ways: it binds the 
substrate less tightly (by ~ 10 2 -fold so) [13] in a groove 
that is significantly wider [7,10]. Lys31 (similar to 
Lys41, its counterpart in RNase A) is critical for catal- 
ysis [13] but most likely requires a substantial (and 
thus far unobserved) change in enzyme conformation 
to participate [7,10]. An additional lysine residue, 
Lys9, also plays a major role in catalysis [13]. Collec- 
tive motion comparable to that of RNase A has not 
been detected, either by essential dynamics simulations 
of the wild-type enzyme [9] or by application of 
dynamic NMR techniques to the <Q1S mutant [27]. 

Our understanding of RNase A and other ribonuc- 
leases has benefited greatly from their characterization 
at atomic or near-atomic resolution [22,28-30]. By con- 
trast, all ONC structure-function studies to date have 
used a significantly lower, 1.7 A resolution structure as 
a reference point [7]. Redressing this, we have deter- 
mined the X-ray crystal structure of ONC at atomic 
resolution. It provides a more precise description of 
the active site region, its geometry and atomic fluctua- 
tions, and allows a refined interpretation of the 
changes that occur upon nucleotide binding [10] and 
mutation [9]. Furthermore, via the combined evalua- 
tion of peptide backbone alignments, anisotropic 
atomic displacement parameters (ADPs) and normal 
mode analysis using an elastic network model, we 
identify collective motion in the ONC structure that 
may have significance for turnover. This provides valu- 
able insight into the interplay between molecular struc- 
ture and collective dynamics during the divergent 
evolution of the pancreatic ribonuclease superfamily. 

Results 

Structural overview 

The atomic resolution structure of ONCSO4 was 
determined from an orthorhombic crystal at 100 K 
(Tables 1 and 2). The previously reported ONCSO4 
structure, determined at room temperature in the same 
crystal form to 1.7 A resolution [Protein Data Bank 
(PDB) entry 10NC] [7] is hereafter referred to as 
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Table 1. Crystallographic statistics. 



Characteristic 


Value 


X-ray wavelength (A) 


0.978 


Resolution range (A) 


32.0-1.1 


Space group 


P2,2 1 2 1 


Unit cell dimensions (A) 


a = 32.23, 




u = oo. id, c = do. yz 


Number of reflections 3 




Measured 


110 620 


Unique 


30 093 


R b 

"symm 


0.059 (0.330) c 


//ct(A 


22.4 (2.8) c 


Completeness (%) b 


85.1 (68.0) c 


Wilson 6-factor 


9.1 


a Data cut-off: / > -3a(/). b B symm = 


XhEill/iW) - </(«>l/S h S i / i (h)], 


where /, is the /th measurement and (1(h)) is the weighted mean of 


all measurements of 1(h). c Values in parentheses refer to the outer- 


most shell (1.12-1.10 A). 




Table 2. Structure refinement statistics. 




Characteristic 


Value 


Number of non-hydrogen atoms 3 




Protein 


842 


Anions 


29 


Water 


149 


Total 


1020 


Number of observations 


29 163 


Number of parameters 


9193 


Number of restraints 


2687 


(Observations + restraints)/parameters 


3.5 


D b 
n work 


0.161 


fifree 0 


0.183 


Mean anisotropy d 


0.57 + 0.14 


Deviation from ideality (rmsd) 




Bond lengths (A) 


0.016 


Bond angles (°) 


1.77 


Mean isotropic equivalent B-factor (A 2 ) 




Protein atoms 


12.1 


Anions 


17.0 


Water molecules 


20.9 


All atoms 


13.4 



a Atoms in alternate conformations are included. b ff wor k = £h 
\F 0 -F 0 \/X h F 0 , where F 0 and F c are the observed and calculated 
structure factor amplitudes of reflection h, respectively. c fi free is 
equal to fl WO rk for a randomly selected 5% of reflections not used 
in the refinement [61]. d Calculated with parvati [69]. 



rtONGS0 4 . The present structure was compared in 
detail with rtONOS0 4 and also with three significant 
low-temperature ONC crystal structures: one in the 
same crystal form [M23L-ONGS0 4 at 1.51 A resolu- 
tion (PDB entry 1YV4)] [9] and two in different crystal 
forms [the ONOd(AUGA) complex at 1.9 A resolu- 



tion (PDB entry 2I5S) and the T89N/E91A-ONG5'- 
AMP complex at 1.65 A resolution (PDB entry 
2GMK )] [10]. Sigma A-weighted electron density maps 
(as supplied by the Uppsala Electron-Density Server) 
[31], as well as the coordinates of all structures, were 
examined. 

The final atomic model fits the experimental data 
well (Table 2). The entire main chain is well defined in 
the 2F 0 -F C electron density map and, of the 95 eligible 
residues, 83 lie in the most favored regions of the 
Ramachandran plot, with the remaining 12 in the 
additional allowed regions. This is similar to the previ- 
ous rtONOS0 4 , M23L-ONGS0 4 and ONGd(AUGA) 
structures but is closer to ideality than the 
T89N/E91A-ONG5'-AMP structure in which 17 resi- 
dues lie in the additional allowed regions. 

The enhanced resolution of the data permitted the 
relaxation of stereochemical restraints during crystallo- 
graphic refinement, providing extra detail throughout 
the structure. A clear illustration of this is provided by 
a segment involving the start of helix H3; here, the 
main chain adopts a strained conformation in which 
the peptide bonds of Ser39, Arg40 and Pro41 exhibit 
significant deviations from planarity (co dihedral 
angles = 160.0, 192.1 and 193.5°, respectively). The 
conformations of the side chains of Arg40 and Glu42 
are well defined, with the guanidino group of Arg40 
coordinating a sulfate ion. These features are not 
apparent in the rtONGS0 4 structure (which was 
refined with a relatively narrow spread of co) and con- 
tribute to a significant shift in the main chain in this 
region, notably a > 1 A change in the position of the 
C a atom of Arg40. 

The data also enabled the resolution of discrete 
alternate conformations for the side chains of six resi- 
dues (Aspl8, Ile51, Ser54, Thr60, Thr83 and Asn84) to 
the extent that their 2F o -F 0 electron density is continu- 
ous when contoured at up to 1.4 a. This contrasts with 
the rtONGS0 4 , M23L-ONGS0 4 and ONC d(AUGA) 
structures, which specify no alternate conformations, 
and with the T89N/E91A-ONG5'-AMP structure, 
where, of the non-mutated residues modeled in alter- 
nate conformations (Asn26, Lys33, Arg40, Glu42 and 
Lys80), only Asn26 has continuous 2F 0 -F C electron 
density, even when the contouring level is reduced to 
1.0 a. 

Active site residues 

The structure provides a refined view of the active site 
as it occurs in the absence of nucleotide bases or sug- 
ars. In accordance with the convention devised for 
RNase A [32], it can be viewed as a series of subsites 
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that bind the successive base and phosphoryl groups of 
the RNA substrate. The scissile phospodiester linkage 
is bound at the Pj subsite, whereas the bases located 
immediately upstream and downstream (5' — > 3') are 
bound at the Bi and B 2 subsites, respectively. 

The majority of active site residues are very well 
defined in the electron density map (Fig. 2A) and 
occupy positions similar to those in the rtONC-S0 4 
structure. Hence, at the Pi subsite, residues Lys9, 
His 10, His97 (the latter in the inactive, B conforma- 
tion) and Phe98 (the counterparts of RNase A residues 
Glnll, His 12, His 119 and Phel20, respectively) coordi- 
nate a sulfate ion (Fig. 2B). The N-terminal pyrogluta- 
mate residue (Peal), which serves to align Lys9 for 
catalysis [13], is now refined in a puckered conforma- 
tion rather than the planar one described by the 
rtONC-S0 4 coordinates. As noted previously [7], 
Lys31 (the counterpart of RNase A Lys41) lies some- 
what further afield. The refinement directs its side 
chain into a more extended conformation than that 
reported for rtONCS0 4 , placing the N c atom 5.22 A 
from the nearest sulfate O atom (cf. 5.37 A in 
rtONCS0 4 ). This leaves it within hydrogen-bonding 
range of the O atom of Lys33 (2.92 A distant) and the 
O 81 atom of Asn34 (3.03 A distant), although it should 
be noted that electron density for Lys31 is lacking 
beyond C s , indicating that the side chain of this 
residue has significant mobility and most likely makes 
no firm hydrogen bonds in the presence of sulfate. 
Reinspection of respective electron density maps indi- 
cates that Lys31 is disordered to a similar extent in 
rtONC-S0 4 but adopts a well-defined, extended con- 
formation in the ONC-d(AUGA) complex. 

There is little positional change among the majority 
of the residues that line the pyrimidine-binding pocket 
(Bi subsite). In addition to HislO and Phe98 noted 
above, this includes Asn34, Thr35, Val99, GlylOO and 
VallOl. The exception is Lys33 (the counterpart of 
RNase A Val43) whose side chain stacks against the 
uracil moiety in the ONCd(AUGA) structure [10]. 
The side chain of this residue is partially disordered 
with weak electron density beyond (Fig. 2A), and 
the refinement suggests that the major side chain rot- 
amer intrudes into the space occupied by the bound 
pyrimidine. In comparison, the rtONCS0 4 structure 
presents the side chain in a more extended conforma- 
tion and the disorder, although detectable, is not as 
pronounced. Upon binding the uracil moiety of 
d(AUGA), the side chain adopts a well-defined, 
extended conformation. 

At the B 2 subsite (a shallow depression lined by 
Thr89, Glu91, Val96 and His97), there is little to dis- 
tinguish the present and the rtONCS0 4 structures. 



Main-chain dynamics 

Comparison with the rtONCS0 4 data reveals that the 
low temperature used during X-ray diffraction brings 
about a 7.7% compaction in unit cell volume, with the 
largest reduction (5.8%) affecting the b axis. This is 
coincident with a compaction of the protein structure, 
as judged by a 2.0% reduction in solvent-excluded vol- 
ume (from 13685 to 13409 A 3 ) and a 1.0% reduction 
in radius of gyration (from 13.736 to 13.603 A). After 
global alignment, the C a traces of the two structures 
differ by only 0.35 A (rmsd). However, difference dis- 
tance analysis [33] of the two structures reveals a wide- 
spread reduction in inter-C a distances (Fig. 3A). 
Furthermore, it indicates that a significant contribution 
to the compaction is made by the movement of lobes 1 
and 2 towards one another. For example, each of the 
strands of the W x P-sheet (B3, B6 and B7, comprising 
residues 55-58, 86-91 and 96-101, respectively) moves 
closer to each of those of V 2 (B0, BI, B4 and B5, com- 
prising residues 11-12, 33-38, 63-70 and 77-84, respec- 
tively). The largest attractive movements are measured 
for the loops and helices that decorate the periphery of 
the lobes, in particular the segments in lobe 2 that con- 
nect B0 to BI (residues 13-32, incorporating helix H2) 
and B4 to B5 (residues 71-75), and the one in lobe 1 
that connects BI to B3 (residues 39-54, incorporating 
helix H3). The conformational change can be visual- 
ized most effectively by aligning the structures on the 
basis of the V! P-sheet residues alone (Figs 4 and 5A). 
This presents an inward movement of the V 2 P-sheet 
and the attached B0-B1 segment that encompasses H2 
and its flanking loops ('segment F). Some of the main 
chain components of the active site cleft are affected. 
For example, the distance between the main-chain N 
atoms of Thr35 and Phe98 (which are located on 
opposite sides of the cleft and have been used to gauge 
its width) [8,9] reduces from 9.17 to 8.86 A. 

Differences in main chain conformation are also 
detectable when the structure is compared with that of 
the low-temperature ONC-d(AUGA) complex. Overall, 
the C a traces of these two structures differ by only 
0.38 A (rmsd) but, again, there are marked deviations 
at the periphery of lobe 2. After a Vi alignment, a 
concerted motion of two of the lobe 2 segments is evi- 
dent (Figs 4 and 5B). One of these is segment I and 
the other comprises the tips of B4 and B5 and the 
intervening loop ('segment IF); the two segments are 
bridged by a disulfide bond (Fig. 1). The segments 
move approximately parallel to the active site cleft in 
the direction of the N-terminus (effectively down the 
nucleotide backbone) but do not bring about further 
compaction of the structure. Indeed, the protein 
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Fig. 2. Active site. (A) Electron density map in the region of the P-i and B, subsites. Shown in stereo are a wireframe representation of the 
final 2F 0 -F C electron density map (contoured at 1 .6 a) and ball-and-stick representations of the protein and the bound sulfate ion. Carbon, nitro- 
gen, oxygen and sulfur atoms are colored gold, blue, red and green, respectively. Spheres denote water molecules and '<Q1' is the N-terminal 
pyroglutamate residue. The majority of residues are very well defined, with the exceptions being Lys31 whose side-shain is disordered beyond 
C s , and Lys33, whose side-chain is poorly defined (occupancy = 0.61) beyond C p . (B) Comparison of sulfate binding by ONC and RNase A. The 
ONC-S0 4 and RNase A-S0 4 (PDB entry 1 KF2 ) [22] structures were superposed on the basis of the C a positions of selected nucleotide-binding 
residues (from ONC, H10, T35 and H97; from RNase A, H12, T45 and H1 1 9). Shown in stereo are ball-and-stick representations of each com- 
plex along with a surface representation of ONC. The color scheme for ONC-S0 4 is the same as for (A), whereas that for RNase A-S0 4 differs 
in that its carbon atoms and sulfate ion are colored pink. Potential hydrogen-bonds are shown as dashed lines (in black for ONC-S0 4 , in pink for 
RNase A-S0 4 ). For those RNase A residues that are observed in dual conformations, the one that most closely matches the conformation of 
the corresponding ONC residue is shown (conformation A of Q1 1 and V43; conformation 8 of K41 , H1 1 9 and the sulfate ion). 



component of the complex has a solvent-excluded vol- 
ume (13430 A 3 ) and radius of gyration (13.636 A) that 
are 0.2% larger than those of the present structure, 



and the difference distance plot indicates an approxi- 
mately even balance of positive and negative C™ dis- 
placements (Fig. 3B). Several residues at the active site 
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Fig. 3. Difference distance analysis of changes in main-chain conformation. Heat maps for difference distance matrices obtained after (A) 
cryocooling (rtONC-S0 4 -> 0NC-S0 4 ) and (B) exchange of a sulfate for a dinucleotide at the active site [0NC-S0 4 -> ONC-d(AUGA)]. Neg- 
ative values (red dots) signify C" pairs that move closer together, whereas positive values (blue dots) signify the converse. 



cleft are affected, notably Lys31, whose C a atom 
moves away from the base at B[ and toward the phos- 
phate at Pi, as well as Lys9 and Thr35, whose C a 
atoms retreat by a smaller margin from the substrate 
analog. The movement reduces the Thr35 N-N Phe98 
distance only slightly to 8.76 A. 



Elements of the aformentioned main-chain dynamics 
are discernible among the ADPs of the present struc- 
ture. The isotropic equivalent ^-factors (B cq ) of the C™ 
atoms are in the range 7.6-15.8 A 2 , with a mean of 
10.3 ± 1.6 A 2 . The majority of the largest B cq values 
are located in lobe 2 with noticeable concentrations in 



Lobe 2 




Lobe 1 



Fig. 4. Main-chain movement in response to cryocooling and dinucleotide binding. C™ traces of ONC extracted from the crystal structures of 
rtONGS0 4 (PDB entry 1QNC ) [7] (gray), 0NC-S0 4 (present study) (blue) and ONC-d(AUGA) (PDB entry 2I5S ) [10] (red), superposed on the 
basis of the V-, (5-sheet residues (black studs). V 2 p-sheet residues are highlighted with red studs, whereas two relatively mobile segments 
(I and II) of lobe 2 are labeled, as are the positions of key substrate-binding residues. To highlight the active site cleft, the protein surface 
and the active site ligand, d(UpGp) (ball-and-stick) from the ONC-d(AUGA) structure are superimposed. The orientation is identical to that of 
Fig. 1. 
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segments I and II (Fig. 5C), confirming that this region 
is a relatively mobile part of the structure. The anisot- 
ropy (A) of the C™ atoms is in the range 0.278-0.856, 
with a mean of 0.610 ± 0.107. Of those C atoms with 



relatively high anisotropy (A below the mean by half 
the standard deviation or more), a high proportion is 
located in segments I and II (Fig. 5C). Here, the prin- 
cipal axes of the C a fluctuations are an excellent match 







Fig. 5. Main-chain dynamics from experiment and theory. In each panel, a C™ trace of ONC is shown, oriented as in Figs 1 and 4 and col- 
ored by an index of mobility (blue represents low; white, intermediate; red, high mobility). (A, B) Information from crystal structure align- 
ment. The rtONC-S0 4 , 0NC-S0 4 and ONC-d(AUGA) crystal structures were superposed on the basis of the V, p-sheet residues, as in Fig. 4. 
(A) Showing rtONC-S0 4 , colored according to the C" displacements resulting from cooling to 100 K (range = 0.02-1.14 A). (B) Showing 
ONC-S0 4 , colored according to the C™ displacements resulting from the exchange of sulfate for d(AUGA) and crystallization in an alternative 
space group (range = 0.04-1.70 A). Arrows represent the C" deformation vectors (up-scaled by an arbitrary factor for clarity). (C) Information 
from the atomic displacement parameters of the ONC-S0 4 crystal structure. The 0NC-S0 4 trace is colored according to the 6 eq values of 
the C" atoms (range = 7.6-15.8 A 2 ). Arrows depict the principal axes of positional fluctuation for those C" atoms whose anisotropy exceeds 
the mean by half the standard deviation or more (i.e. A < 0.5565), scaled in proportion to the respective eigenvalue. (D-F) The three lowest 
frequency nontrivial modes of motion (numbers 7-9, respectively) predicted by ANM. Traces of 0NC-S0 4 are colored according to the 
relative mobility at each position; arrows show the direction and magnitude of motion during the compaction phase of each mode (arbitrary 
scaling). 
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for the axes of the displacements deduced from struc- 
tural alignment with the ONOd(AUGA) complex 
(Fig. 5B), indicating that, after the binding of sulfate, 
these segments have an intrinsic propensity to undergo 
this motion. 

Intrinsic collective motion was further investigated 
using the anisotropic network model (ANM) with nor- 
mal mode analysis. The three lowest-frequency non- 
trivial modes (numbers 7-9) are found to describe 
variations on the compaction/expansion ('breathing') 
theme (Fig. 5D-F). In each case, there is a counter- 
motion of the two lobes, with lobe 2 undergoing the 
larger movement. Modes 7 and 9 also feature a propel- 
ler-like twisting of lobe 2, such that segments I and II 
move in opposite directions relative to the active site 
cleft. When viewed during the equivalent phase of 
motion (e.g. compaction), it is evident that the twist 
directions in these two modes are opposed. Mode 8, 
on the other hand, describes a pincer-like counter- 
motion of the two lobes in which segments I and II 
move in unison. This is reminiscent of the breathing 
motion identified in RNase A via crystallographic 
analyses [8,20,23,24] and molecular dynamics simula- 
tions [34-36]. 

The compaction phases of modes 8 and 9 have much 
in common with the response of crystalline ONC S0 4 
to cryocooling (Fig. 5A), although they both overesti- 
mate the movement of segment II. A plausible expla- 
nation for this lies in the intimate involvement of this 
segment in crystal packing interactions that constrain 
its movement. The additional movement of lobe 2 that 
accompanies nucleotide binding (Fig. 5B) may have its 
origins in multiple modes. For example, there are simi- 
lar trends for segment II in the compaction phase of 
mode 9, for segment I in the compaction phase 
of mode 7, and both segments in the expansion phase 
of mode 8. Interestingly, when comparable ANM com- 
putations are applied to the structure of RNase A 
(PDB entry 7RSA ) [28], all three types of breathing 
motion are detected (data not shown). 

Comparison with M23L-OIMC S0 4 

The structural details of superactive ONC mutants 
may offer valuable insight into the low ribonucleolytic 
activity of the wild-type enzyme. One such mutant, 
M23L-ONC, exhibits a five-fold elevation in RNA 
cleavage activity, a significant reduction in thermosta- 
bility and an increase in susceptibility to proteolysis 
[37-39]. The C a trace of M23L-ONCS0 4 as deter- 
mined crystallographically at 100 K (PDB entry 1YV4 ) 
[9] bears a very strong resemblance to that of the pres- 
ent structure; the overall C a deviation is 0.15 A (rmsd), 



with the greatest deviation (0.37 A) occuring at the 
mutation site. Using the rtONCS0 4 structure as a ref- 
erence point, Merlino et al. [9] identified, in M23L- 
ONCS0 4 , a difference in the conformation of Lys31 
and a change in the hydrogen-bonding partner of its 
N 5 atom (from Lys33 O in the wild-type to Asn34 O 81 
in the M23L mutant). They suggested that this placed 
it in a more favorable position to activate the enzy- 
matic reaction [9]. The present structural analyses 
refine our view of these features: the disorder that we 
have observed in the Lys31 side-chain indicates that it 
is unlikely that a Lys31 >r— O Lys33 bond (or indeed 
any hydrogen bond involving Lys31 N^) exists in 
ONOSO4 at room temperature or at 100 K. An 
important effect of the mutation, therefore, is the 
rigidification of the Lys31 side-chain. If viewed as a 
reduction in ground-state entropy, this represents one 
possible means by which catalysis is enhanced in this 
mutant. An additional means stems from the apparent 
overall increase in conformational flexibility (as indi- 
cated by the reductions in thermostability and resis- 
tance to proteolysis), which may favor the significant 
deformation necessary to attain the transition state. 

Elsewhere across the substrate-binding site, there is 
little difference between the present structure and that 
of M23L-ONCS0 4 . This extends to the disorder in the 
side-chain of Lys33. 

Discussion 

It has become increasingly clear that the primary 
sequence of a protein encodes not only its equilibrium 
3D structure, but also its intrinsic collective motion 
[40-42]. Multidisciplinary studies of model enzyme sys- 
tems such as adenylate kinase have shown that these 
large-scale fluctuations in conformation can be coupled 
to catalysis [43]. It is likely that this is the case for 
RNase A, in which the rate constant for the predomi- 
nant collective motion (1700 s -1 ) is the same as both 
k CSLt and the rate constant for product dissociation 
[24-26]. Because patterns of backbone flexibility (in 
particular, the largest-scale collective modes) tend to 
be conserved within structurally homologous protein 
families [40-42], a proposed general conformational 
reaction scheme for catalysis by members of the pan- 
creatic ribonuclease superfamily is presented in Fig. 6. 
Consistent with the behavior of RNase A, conforma- 
tional minima (E n ) that are intrinsic to the enzyme are 
mapped across the reaction sequence. In the absence 
of substrate, the enzyme samples a wide range of con- 
formations (including those on the catalytic pathway), 
with the equilibria strongly favoring E\ (an open 
conformation that corresponds to the global energy 
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Fig. 6. General scheme for the coupling of conformational change 
and catalysis in pancreatic ribonuclease superfamily members. E, is 
the global conformational minimum and E 2 , E 3 ... E n are local con- 
formational minima along the catalytic pathway. After binding of 
the substrate (S) and attainment of the transition state (EjS*), the 
reaction yields two products (P). In the case of a dinucleotide sub- 
strate (i.e. NpN'), the products are N 2',3'-cyclic phosphate and N'. 
Release of the two products may involve multiple steps, as 
denoted by the dashed arrows. 

minimum). The substrate is therefore, most likely to 
bind to E x but, because the E x — > E 2 conformational 
change is only a subtle closure, it could conceivably 
also bind to E 2 or any of the conformational micro- 
states inbetween. Substrate binding shifts the equilib- 
rium to favor E 2 , from which the enzyme undergoes a 
further conformational change (to E 3 ) to optimize 
complementarity to the transition state, most likely 
involving a more extreme closure [44]. Subsequent con- 
formational changes effect release of the two products 
and complete the catalytic cycle. It follows that, when 
trying to understand the distinctive enzymatic proper- 
ties of ONC in relation to those of RNase A, an 
attempt must be made to compare both the local 
active site structure and the relevant collective dynam- 
ics of each enzyme. 

Substrate-binding residues 

A deficit in nucleotide-binding affinity accounts for a 
large proportion of the low ribonucleolytic activity of 
ONC relative to that of RNase A [13]. Previous crys- 
tallographic work on the binding of (deoxy)dinucleo- 
tides to the active site of each enzyme has indicated 
that enzyme-nucleotide interactions are not well opti- 
mized when these compounds initially bind to ONC 
[10,15]. At the Pi subsite, there is one less protein 
atom within hydrogen-bonding range of the phospho- 
diester linkage than at the same subsite of RNase A 
(in ONC: Lys9 N<\ His 10 N e2 and Phe98 N; in 
RNase A: Glnll N e2 , Hisl2 N e2 , Hisll9N 51 and 
Phel20 N), whereas the atom of Lys31 is 1.3 A fur- 
ther away from the nearest phosphoryl oxygen than is 
its counterpart, Lys41 N^, in RNase A (in ONC: 
5.0 A; in RNase A, 3.7 A). Indeed, the side-chain of 
Lys31 in ONC does not contribute to ground-state 
nucleotide binding [13], unlike the side-chain of Lys41 



in RNase A [45]. The findings of the present study 
imply that, despite the lack of a direct or water-medi- 
ated interaction between Lys31 and the nucleotide, the 
binding process brings about a rigidification of this 
residue. A comparable reduction in entropy is 
observed for the side-chain of Lys41 when nucleotide- 
free RNase A crystal structures (e.g. PDB entries 
7RSA [28] and 1KF2 [22]) are compared with the 
RNase A d(CpA) complex (PDB entry 1RPG [15]). By 
similar deduction, a rigidification of Lys33 at the B ; 
subsite of ONC also accompanies nucleotide binding. 
Although the equivalent residue in RNase A (Val43) 
shows the same tendency, its disorder is lower, being 
limited to discrete alternate conformations in the 
absence of nucleotides. Indeed, Val43 has substantially 
fewer degrees of freedom than does Lys33; thus, it is 
possible that an increased entropic penalty at this posi- 
tion contributes to the reduced nucleotide-binding 
affinity of ONC. 

Crystallographically-observed conformational 
variations 

The reaction scheme in Fig. 6 demands that the 
enzyme possesses the requisite degree of large-scale 
conformational flexibility. In RNase A, this has been 
established experimentally via analysis of multiple 
crystal structures [23] and us/ms-timescale solution 
dynamics [26]. In comparison, the number of available 
ONC crystal structures is small and no solution struc- 
ture/dynamics of wild-type ONC have been reported. 
Our analyses of the present structure have, for the first 
time, identified appreciable collective motion in ONC. 
In all the analyses, the deduced motion has some 
similarity to that identified in RNase A in that it is 
greatest in the loops and helices that decorate the 
P-structure. Most of these loops are shorter in ONC 
than in RNase A (amounting to 20 fewer residues 
overall) and the locations of the longest loops appear 
to bias the motion toward one side of the protein 
(lobe 2). The magnitude of the deformation that occurs 
in response to cryocooling is modest but comparable 
to that measured for RNase A [19], whereas the effect 
of dinucleotide binding on these loops (in conjunction 
with crystallization in an alternative space group) is 
far more pronounced, and indeed more than when 
d(CpA) binds to RNase A [15]. 

The impact of the observed collective motion on the 
width of the active site cleft is of relevance to its poten- 
tial role in catalysis. In the E x conformation, the cleft 
is most likely significantly wider than that in RNase A 
[9]. To clarify how much wider, it would be helpful to 
compare the unliganded conformations of ONC and 
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RNase A under closely matched conditions. Although 
suitable crystal and solution structures of RNase A are 
available [28,46,47], none of the ONC crystal structures 
reported to date (including the present structure) are 
unliganded and none of the ONC NMR structures are 
of the wild-type enzyme. In the absence of these data, 
a useful guide comes from comparison of the room- 
temperature, sulfate-liganded structures at pH 5.2-5.5. 
Using, as metrics, the Thr35 N-N Phe98 distance 
in ONC and the Thr45 N-N Phel20 distance in 
RNase A [8,9], the width of the cleft in ONC (9.17 A; 
PDB entry 1QNC ) [7] is ~ 0.8 A more than that in 
RNase A (8.36 A; PDB entry 1KF2) [22]. It must be 
noted that even the binding of sulfate (a minimal sub- 
strate analogue) has a measurable effect on the global 
conformation of RNase A, reducing the cleft width of 
the unliganded enzyme by ~ 0.2 A (8.55 A; PDB entry 
7RSA ) [28]. Unliganded ONC, therefore, may present 
a substrate-binding cleft with a width in excess of 
9.17 A. 

The narrowing of the active site cleft caused by cool- 
ing ONGSO4 crystals from room temperature to 
100 K (0.31 A) is comparable to that observed after 
cooling the P2, crystal form of RNase A from 300 to 
98 K (0.28 A; PDB entries 8 RAT and 1RAT , respec- 
tively) [19]. In addition, the marginal effect on cleft 
width observed on exchanging sulfate for d(UpGp) at 
the active site of ONC at constant temperature (a 
decrease of 0.10 A) matches exactly that observed with 
sulfate- and d(CpA)-liganded RNase A (PDB entries 
1KF2 and 1RPG , respectively) [15,22]. Hence, the pos- 
sibility that the width of the active site cleft of ONC is 
modulated by an RNase A-like breathing motion 
remains open. A stronger stimulus of conformational 
change may allow a more extreme conformational of 
ONC to be trapped crystallographically or in solution. 
In the case of RNase A, the effects of the two environ- 
mental factors evaluated in the present work (tempera- 
ture and dinucleotide binding) are surpassed by those 
resulting from the binding of pyrimidine nucleotide 
phosphates such as 2'-UMP [48] and 2'-CMP [8,15], 
which can be considered as product analogs. An even 
greater conformational change is likely to result from 
the binding of a transition state analog [44]. Indeed, a 
structural model for the ES* complex as it occurs in 
RNase A or ONC would constitute a major advance- 
ment in our understanding of the structural and 
dynamic requirments for catalysis by these enzymes. 
The crystallographically-observed binding mode of uri- 
dine vanadate, a marginal transition state analog [45], 
was not consistent with that of a transition state [49], 
so it is important that other transition state analogue 
candidates are evaluated. 



Although one might hope to trap the 'landmark' 
conformations depicted in Fig. 6 via crystallization, it 
must be noted that each is surrounded by a large num- 
ber of conformational microstates [50]. Because of the 
geometrical constraints of crystal packing, some of 
these may be much more conducive to lattice forma- 
tion/preservation, leading to their prevalance over the 
landmark conformation. The influence of the lattice 
can be strong, perhaps most clearly illustrated by the 
C2 crystal form of RNase A in which one of the two 
protomers in the asymmetric unit (mol B) is con- 
strained in a more open conformation than the other, 
irrespective of the ligand that is bound. For example, 
the active site cleft of mol B when bound to 2'-UMP is 
0.65 A wider than that of mol A in the same state, 
and is 0.53 A wider than that of mol A in unliganded 
form (PDB entries 1O0M [47] and 2G8Q [48], respec- 
tively). Hence, it would be useful to study multiple 
crystal forms of ONC before drawing firm conclusions 
on its conformational flexibility. 

Intrinsic collective motion 

In view of the present results and the tendency for the 
conservation of collective modes of backbone motion 
within structurally homologous protein families [40- 
42], it is most likely that the differences in sequence 
between ONC and RNase A have altered the relative 
populations of different conformational states and/or 
the kinetics of their interconversion [50]. NMR relaxa- 
tion experiments on <QlS-ONC suggest that the 
wild-type enzyme is conformationally more uniform 
than is RNase A [27], so the conformational landscape 
sampled by unliganded ONC may favor the open con- 
formation (Ei) to such an extent that the other confor- 
mational states are invisible to the technique. This 
need not be a large change from the equilibrium posi- 
tion measured for RNase A: it was calculated that a 
3% upward shift in the relative population of E\ from 
its wild-type value of 95% would be sufficient to 
remove the NMR relaxation dispersion signal [51]. In 
addition, because the binding of a substrate or product 
analog to RNase A shifts the equilibrium to favor the 
respective closed conformation [26], NMR relaxation 
experiments conducted with ONC in the presence of 
such analogs may be more likely to reveal alternative 
conformational states [52]. 

Despite their simplicity, coarse-grained elastic net- 
work models such as the ANM have proven to be 
remarkably robust in their prediction of global protein 
dynamics [53]. The ANM analysis in the present study 
indicates that the largest-scale (lowest-frequency) 
collective modes of motion are conserved between 
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RNase A and ONC. This includes three compaction/ 
expansion variations, one of which resembles the 
familiar pincer-like breathing motion of RNase A. 
Although this breathing motion appears to predomi- 
nate overwhelmingly in RNase A (as portrayed across 
numerous crystallographic 'snapshots' [23] and in the 
anisotropy of its high-resolution structure [54]), the 
findings of the present study suggest that this may not 
be the case for ONC. It is possible that, as a result of 
the more open E\ conformation of ONC, a more com- 
plex conformational change is required to effect cataly- 
sis. The ANM results are in contrast to the essential 
dynamics extracted from all-atom molecular dynamics 
simulations of 3-10 ns in duration, which highlighted a 
breathing motion in RNase A [23,36] but not in ONC 
[9]. Although atomistic molecular dynamics is consid- 
ered to be the 'gold standard' for modeling protein 
fluctuations, the lowest-frequency modes of motion 
can be significantly under-represented if the simulation 
time is less than several hundred nanoseconds [55]. 
Therefore, it may be that longer molecular dynamics 
simulations are required to detect large-scale motion in 
ONC. 

ONC now joins a select group of ribonucleases 
(comprising RNase A [22], RNase Sa [29] and eosino- 
phil-derived neurotoxin [30]) for which an atomic reso- 
lution structure has been determined. The structure 
reported in the present study has enabled the re-inter- 
pretation of existing work and will serve as a valuable 
reference point for future structure-function studies on 
the protein. 

Note added in proof 

After the submission of the present manuscript, the 
authors became aware of two unpublished ONC struc- 
tures that have been deposited in the PDB. These are 
both low-temperature, sulfate-liganded structures in 
the same crystal form as the present structure (PDB 
entries 3HG6 and 3PHN at 1.70 and 1.46 A resolu- 
tions, respectively). They exhibit modest conforma- 
tional differences from the present structure, the 
nature of which are in accordance with the conclusions 
drawn above. 

Experimental procedures 

Crystallization and data collection 

Authentic <Ql-ONC was purified from R. pipiens oocytes 
by a modification [56] of the original method [1]. Crystals 
were grown at 16 °C by the hanging-drop/vapor-diffusion 



method. In a refinement of previously established condi- 
tions [7], 2 u.L of protein solution (10 mgmL" 1 in water) 
was mixed with an equal volume of reservoir solution 
[2.2 M (NH 4 ) 2 S0 4 , 0.25% (w/v) polyethylene glycol) 4000, 
0.01 M sodium acetate, pH 4.5], producing a rod-like crystal 
within 1 month. The crystal was flash-cooled to 100 K 
(using 25% glycerol) and diffraction data were collected on 
station 14.2 of the Synchrotron Radiation Source (Dares- 
bury, UK). Images were processed using the hkl-2000 suite 
[57], as summarized in Table 1. The crystal was isomor- 
phous with those grown under the parental conditions [7], 
containing one protein molecule per asymmetric unit. 

Structure determination and refinement 

Structure factor amplitudes were calculated using the trun- 
cate procedure [58,59] and initial phases were obtained 
using the molecular replacement routines of amore [60], 
employing, as a search model, the ONC coordinates from 
PDB entry 1QNC [7]. All subsequent refinement adopted a 
maximum-likelihood target function and used a randomly- 
selected 5% of reflections for cross-validation [61]. Initial 
rounds of refinement were carried out with REFMAC5 [62] 
and were interspersed with electron density map viewing 
and manual model adjustments using coot [63]. During this 
time, isotropic ADPs were used. After inclusion of anions 
and water molecules, it was verified by sfcheck [64] that 
the diffraction data had no significant angular bias within 
the reciprocal lattice, anisotropic ADPs were adopted and 
then refinement proceeded with phenix [65]. The ADP 
X-ray weighting term wxujscale was maintained at the 
default value of 1.0, whereas the stereochemical X-ray 
weighting term wxc_scale was increased from 0.5 to 1.5. 
For residues whose side-chains were observed in discrete 
alternate conformations, the occupancy of each conforma- 
tion was refined (sum of occupancies = 1). The occupan- 
cies of three side-chains (Lys33, Lys78, Lys80) and three 
sulfate ions with single, weakly-defined conformations were 
also refined (occupancy < 1). In the final stages, hydrogen 
atoms were added in riding model positions. Structure vali- 
dation was conducted with procheck [66], what^check 
[67], molprobity [68] and parvati [69]. Detailed refinement 
statistics are given in Table 2. 

Structural analysis 

Electron density maps for previously published protein 
structures were obtained from the Uppsala Electron-Den- 
sity Server [31]. Structural alignments were performed with 
vmd [70], as were mass-weighted radius of gyration calcula- 
tions. Solvent-excluded volume was calculated using the 
Voronoi polyhedra approach [71] with a probe radius of 
1.4 A as implemented in vadar [72]. Before calculation of 
the solvent-excluded volume and radius of gyration, 
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secondary side-chain conformations, hydrogen atoms, 
solutes and water molecules were removed from coordinate 
files. Difference distance matrix plots were generated using 
ddmp (Center for Structural Biology, Yale University, New 
Haven, CT, USA). 

Anisotropic ADPs were analyzed using the prody pack- 
age [73]. Briefly, for each C* atom, the ADP matrix was di- 
agonalized to obtain the three principal axes (eigenvectors) 
and their mean-square fluctuations (eigenvalues). At each 
position, the anisotropy, A was defined as the ratio of the 
smallest to the largest eigenvalue [69]. Intrinsic collective 
motion was investigated using the ANM [74] as imple- 
mented in prody. ANM calculations involved only C™ 
atoms and employed a uniform spring constant with a cut- 
off distance of 15 A for interacting atoms. The six slowest 
nontrivial normal modes were examined. Molecular graph- 
ics were produced using vmd, nmwiz [73] and CCP4MG [75]. 
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